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We present new experimental constraints on the elastic, spin-dependent WIMP-nucleon cross 
section using recent data from the XENON100 experiment, operated in the Laboratori Nazionali 
del Gran Sasso in Italy. An analysis of 224.6 live days x 34 kg of exposure acquired during 2011 and 
2012 revealed no excess signal due to axial-vector WIMP interactions with 129 Xe and 131 Xe nuclei. 
This leads to the most stringent upper limits on WIMP-neutron cross sections for WIMP masses 
above 6GeV/c 2 , with a minimum cross section of 3.5xl0 _40 cm 2 at a WIMP mass of 45GeV/c 2 , 
at 90% confidence level. 
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XENON100 was built to search for hypothetical, 
weakly interacting massive particles (WIMPs), which 
could explain the non-baryonic, cold dark matter in 
our Universe [1]. Independently of astrophysical and 
cosmological observations, WIMPs are a consequence 
of many extensions of the Standard Model of particle 
physics, as new, stable or long-lived neutral particles. 
The WIMP dark matter hypothesis is testable by ex- 
periment, the most compelling avenue is to directly ob- 
serve WIMPs scattering off atomic nuclei in ultra-low 
background terrestrial detectors [2, 3]. XENON100 is a 
double-phase xenon time projection chamber operated 
at the Laboratori Nazionali del Gran Sasso (LNGS) in 
Italy. A total of 178 low- radioactivity, UV-sensitive pho- 
tomultiplier tubes detect the prompt (SI) and propor- 



tional (S2) light signals induced by particles interacting 
in the sensitive volume, containing 62 kg of ultra-pure 
liquid xenon. The background level in the energy re- 
gion of interest for dark matter searches (<50keV nr ) 
is 5.3 xlO -3 events kg^ 1 d" 1 keV^ 1 , before discrimination 
of electronic and nuclear recoils based on their S2/S1- 
ratio [4, 5]. The instrument is described in [6], the anal- 
ysis procedure is detailed in [7]. 

WIMPs in the halo of our Galaxy are expected to 
be highly non-relativistic [8] and their interactions with 
nuclei can be characterized in terms of scalar (or spin- 
independent, SI) and axial- vector (or spin-dependent, 
SD) couplings [1, 2]. In the case of SI interactions, the 
leading contribution of the scattering is coherent across 
the nucleus, and roughly scales with A 2 , where A is 



2 



the number of nucleons. Our SI result was presented 
in [4] and excludes a WIMP-nucleon cross section above 
2 x 1CT 45 cm 2 at a WIMP mass of 55 GeV/c 2 at 90% 
confidence level. Here we use the same data set, with 
an exposure of 224.6 live days, a fiducial mass of 34 kg, 
identical event selection cuts, acceptances, relative scin- 
tillation efficiency and background model to derive limits 
on spin-dependent interactions. 

If the WIMP is a spin- 1/2 or a spin-1 field, the con- 
tributions to the WIMP-nucleus scattering cross section 
arise from couplings of the WIMP field to the quark axial 
current. In the case of the lightest neutralino in super- 
symmetric models for instance, scattering occurs through 
the exchange of Z bosons or squarks [1]. To predict ac- 
tual rates, these fundamental interactions are first trans- 
lated into interactions with nucleons by evaluating the 
matrix element of the quark axial-vector current in a 
nucleon. Finally, the spin components of the nucleons 
must be added coherently using nuclear wave functions 
to yield the matrix element for the SD WIMP-nucleus 
cross section as a function of momentum transfer. The 
SD differential WIMP-nucleus cross section as a function 
of momentum transfer q can be written as [9]: 



d<J S p(q) 
dq 2 
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where Gf is the Fermi constant, v is the WIMP speed 
relative to the target, J is the total angular momentum of 
the nucleus and Sa is the axial-vector structure function. 
In the limit of zero momentum transfer [32] the structure 
function reduces to the form [101: 



S A (0) = 



(2J+ 1)(J+ 1) 
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where (S p , n ) = (J\S Ptn \J) are the expectation values of 
the total proton and neutron spin operators in the nu- 
cleus, and the effective WIMP couplings to protons and 
neutrons are defined in terms of the isoscalar ao = a p +a n 
and isovector a\ = a p — a n couplings. 

WIMPs will thus couple to the total angular momen- 
tum of a nucleus and only nuclei with an odd number 
of protons or /and neutrons will yield a significant sensi- 
tivity to this channel. Natural xenon contains two non- 
zero spin isotopes, 129 Xe (spin-1/2) and 131 Xe (spin-3/2), 
with an abundance of 26.4% and 21.2%, respectively. In 
XENON100, the isotopic abundances of 129 Xe and 131 Xe 
are changed to 26.2% and 21.8%, respectively. 

To compare results from different target materials, a 
common practice is to report the cross section for the 
interaction with a single nucleon (a pi a n ) [11—13]. As- 
suming that WIMPs couple predominantly to protons 
(a n = 0) or neutrons (a p = 0), the WIMP-nucleon cross 
section becomes: 
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where <tsd is the total WIMP-nucleus cross section, [ia 
and // Pj71 are the WIMP-nucleus and WIMP-nucleon re- 
duced masses, respectively. 

Calculations of the structure functions Sa(q) are tra- 
ditionally based on the nuclear shell model, but dif- 
fer in the effective nucleon-nucleon interactions and in 
the valence space and truncation used for the computa- 
tion. For xenon as a WIMP target material, we consider 
three large-scale shell-model calculations: by Ressell and 
Dean [14] with the Bonn-A [15] two-nucleon potential, by 
Toivanen et al. [16], using the CD-Bonn potential [15], 
and the recent results by Menendez et al. [17], using 
state-of-the-art valence shell interactions [18, 19] and less 
severe truncations of the valence space. Menendez et al. 
[17] also use for the first time chiral effective field theory 
(EFT) currents [20] to determine the couplings of WIMPs 
to nucleons. The currents for spin-dependent scattering 
are derived at the one-body level and the leading long- 
range two-nucleon currents are included, resulting in a 
reduction of the isovector part of the one-body axial- 
vector WIMP currents [17]. The resulting chiral EFT 
currents are then used to calculate the structure func- 
tions for the WIMP-xenon scattering. Theoretical errors 
due to nuclear uncertainties can be provided when chiral 
two-body currents are included [17]. We show their ef- 
fect in Figure 3. The shell-model calculations are based 
on the largest many-body spaces accessible with nuclear 
interactions, also used to calculate double-beta decay ma- 
trix elements for nuclei up to 136 Xe and to study nuclear 
structures [18, 19, 21]. 

The new calculations by Menendez et al. [17] yield a 
far superior agreement between calculated and measured 
spectra of the 129 Xe and 131 Xe nuclei, both in energy and 
in the ordering of the nuclear levels, compared to older 
[16] results. The values for (S Pt n) ar e close to those of 
Ressell and Dean [14], but quite different from the re- 
sults of Toivanen et al. [16], as summarized in Table I. 
We thus use the Menendez et al. [17] structure functions 
for our benchmark upper limits on WIMP-neutron and 
WIMP-proton cross sections. We also provide a compar- 
ison to the limits obtained when using the calculations 
by Ressell and Dean [14] and Toivanen et al. [16]. In all 
cases, |(5„)| S> | (S^) | , as expected for the two xenon nu- 
clei with an odd number of neutrons and an even number 
of protons. 

Figure 1 shows the structure functions Sa(q) obtained 
from the three calculations for pure neutron and pure 
proton couplings as a function of nuclear recoil energy. 
For the neutron coupling case, for which xenon has the 
highest sensitivity, the functions are rather similar. For 
the proton coupling case, the structure function by Toiva- 
nen et al. [16] differs significantly from the other two re- 
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TABLE I: Parameters of the xenon isotopes used in this analysis: nuclear total angular momentum and parity of the ground 
state, J F , predicted expectation values of the total proton and neutron spin operators in the nucleus (S n ,p) by the Ressell and 
Dean (Bonn A potential) [14], Toivanen et al. (Bonn CD potential) [16] and Menendez et al. (state-of-the art valence shell 
interactions) [17] calculations. 



Ressell and Dean [14] 


Toivanen et al. [16] 


Menendez et al. [17] 


Nucleus 


J p 


(Sn) (S p ) 


(Sn) (Sp) 


(Sn) (S p ) 


129 Xe 


\2) g. S . 


0.359 0.028 


0.273 -0.0019 


0.329 0.010 


131 Xe 


\2)g.s. 


-0.227 -0.009 


-0.125 -0.00069 


-0.272 -0.009 




Recoil 



FIG. 1: Structure functions for 129 Xe (top) and 131 Xe (bot- 
tom) for the case of neutron (plain) and proton (dashed) cou- 
plings, as a function of recoil energy using the calculations of 
Ressell and Dean [14], Toivanen et al. [16] and Menendez et 
al. [17]. The difference is most significant in the case of the 
proton coupling for the Toivanen et al. results. 

suits. Table I summarizes the expectation values of the 
total proton and neutron spin operators in the nucleus for 
129 Xe and 131 Xe in the zero momentum transfer limit. 

Constraints on the spin-dependent WIMP-nucleon 
cross sections are calculated using the Profile Likelihood 
approach described in [31]. Systematic uncertainties in 
the energy scale and in the background expectation are 
taken into account when constructing the Profile Like- 
lihood model and are reflected in the actual limit. It 
is given at 90% C.L. after taking into account statisti- 
cal downward fluctuations in the background. We as- 
sume that the dark matter is distributed in an isothermal 
halo with a truncated Maxwellian velocity distribution 
with a local circular speed of v c — 220km/s, galactic 
escape velocity v csc — 544km/s and a local density of 
p = 0.3GeVcm- 3 [8]. 

The resulting upper limits from XENON100, along 
with results from other experiments, are shown in Fig- 
ure 2 for neutron couplings (top panel) and proton cou- 
plings (lower panel). The la (2cr) uncertainty on the 




FIG. 2: XENON100 90% C.L. upper limits on the WIMP 
SD cross section on neutrons (top) and protons (bottom) us- 
ing Menendez et al. [17]. The la (2a) uncertainty on the 
expected sensitivity of this run is show as a green (yellow) 
band. Also shown are results from XENON10 [22] (using 
Ressel and Dean [14]), CDMS [23, 24], ZEPLIN-III [25] (us- 
ing Toivanen et al. [16]), PICASSO [26] , COUPP [27], SIM- 
PLE [28], KIMS [29], IceCube [30] in the hard {W+W , t+t' 
for WIMP masses <80.4 GeV/c 2 ) and soft (66) annihilation 
channels. 
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-Toivanen et al. 
Mendendez et al. 
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interpretation of the results in terms of SD pure-proton 
cross section strongly depends on the used nuclear model. 
However, regardless of the nuclear model, we obtain the 
most stringent limits to date on spin-dependent WIMP- 
neutron couplings for WIMP masses above 6GeV/c 2 at 
90% C.L. 
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merical form. We thank J. Suhonen for providing us the 
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FIG. 3: The ratio of upper limits calculated with the Ressell 
and Dean [14] and the Toivanen et al. [16] results for the 
structure functions to the ones obtained using Menendez et 
al. [17] for the case of neutron (top) and proton (bottom) 
couplings, along with the theoretical uncertainty band due to 
chiral two-body currents [17]. 



sensitivity of this run, namely the expected limit in ab- 
sence of a signal above the background, is shown as 
a green (yellow) band in Figure 2. The impact on 
these limits when using the Toivanen et al. and the 
Ressell and Dean calculations are shown in Figure 3. 
XENON100 provides the most stringent limits for pure 
neutron couplings for WIMP masses above 6GeV/c 2 , ex- 
cluding previously unexplored regions in the allowed pa- 
rameter space. The minimum WIMP-neutron cross sec- 
tion is 3.5xlCr 40 cm 2 at a WIMP mass of 45GeV/c 2 , 
using Menendez et al. [17]. It changes to 2.5xl0~ 40 cm 2 
and 4.5x 10 _40 cm 2 when using Ressell and Dean and 
Toivanen et al., respectively. The sensitivity to proton 
couplings (Figure 2, bottom panel) is much weaker be- 
cause, as detailed above, both 129 Xe and 131 Xe have an 
unpaired neutron but an even number of protons, thus 
\(Sp}\ <C \(S n }\ (see table I). Upper limits from other 
direct and indirect detection experiments are shown for 
comparison. 

In conclusion, we have analyzed data from 224.6 live 
days x 34 kg exposure acquired by XENON100 during 13 
months of operation in 2011/2012 for SD WIMP interac- 
tions. We saw no evidence for a dark matter signal and 
have obtained new experimental upper limits on the spin- 
dependent WIMP-nucleon cross section. For our limits, 
we use the new calculations by Menendez et al. [17], 
where the WIMP couplings to nucleons are derived us- 
ing chiral EFT currents and which yield a good agree- 
ment between the calculated and measured energy spec- 
tra of the 129 Xe and 131 Xe nuclei. We note that the 
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